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discharge 


G.I.  Sukhinin.  A.V.  Fedoseev 

Institute  of  Thermophysics,  Novosibirsk,  630090,  Russia 

A spherical  glow  discharge  with  a central  point-like  anode  has  been  studied  in  self-consistent 
kinetic  approach.  The  model  includes  non-local  Boltzmann  equation  coupled  with  non-stationary 
ion  balance  equation , and  Poisson  equation  for  radial  electric  field.  For  low-pressure  argon,  the 
model  permits  to  obtain  moving  striations  in  a spherical  glow  discharge . 


1.  Introduction 

Spherical  striations  were  observed  in  low-pressure  glow 
discharge  around  point-like  anode  [1].  It  was  found  that 
for  different  discharge  conditions  the  radii  of  striations 
in  a positive  column  of  discharge  obey  relationship 
fk+i/rk  * const  (geometric  series)  or  rk+]  ~ rk  « const 
(arithmetical  progression).  The  reason  for  such  different 
behavior  of  spherical  striations  is  not  clear  yet.  In  [2], 
spherical  striations  were  modeled  in  a given  peak- like 
electric  field.  However,  the  solution  of  [2]  was  not  self- 
consistent. 

In  the  paper,  this  problem  is  studded  on  the  basis  of 
non-local  Boltzmann  equation  for  electron  distribution 
function  coupled  with  non-stationary  balance  equation 
for  ions,  and  Poisson  equation  for  self-consistent 
electric  field. 


2.  Model 

The  model  includes: 

1.  The  non-local  Boltzmann  equation  in  two-term 
approximation  for  isotropic  part  of  electron  energy 
distribution  function  (EEDF)  written  in  “total  energy- 
space  coordinate  (radius  from  point-like  anode)” 
representation:. 
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where  Ng  is  gas  density,  U is  kinetic  energy  of  electron, 
8=U-e0<p(r)  is  total  electron  energy,  <p(r)  is  radial 
dependence  of  spherical  discharge  potential  in  the 
positive  column  (PC),  £(r)=-3(p/9r  is  electric  field  in  the 
positive  column  of  a spherical  discharge.  Eq.  (1)  is 
parabolic  equations  with  source  terms  representing 
elastic,  inelastic  and  ionizing  electron-atom  collisions. 

2.  The  non-stationary  continuity  equation  in  drift 
approximation  for  ions 
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Vj(r)  is  the  direct  ionization  frequency,  \i\  is  ion  mobility 
coefficient. 


3.  Poisson  equation  for  self-consistent  electric  field: 
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The  solution  of  the  system  (1-3)  was  looked  for  in  the 
region  from  the  edge  (r=R)  of  the  positive  column  to  the 
central  anode  with  appropriate  boundary  conditions. 


The  iterative  numerical  procedure  for  the  solution  of  the 
striation  problem  is  used.  At  the  first  step,  in  order  to 
obtain  the  first  approximation  of  distribution  function 
//(r,U),  the  parabolic  equation  (1)  is  solved  in  some 
proposed  electric  field  E°(r).  From  the  //(r,U),  the 
electron  density  distribution  ne°(r)  and  the  ionization 
frequency  Vi°(r)  are  obtained.  With  the  help  of  Eqs.(2- 
3),  new  approximation  of  electric  field  E(r)  is  found  and 
the  procedure  is  repeated  until  successive  iterations  of 
electric  field  are  converged.  Usually  ten  iterations  are 
enough  for  the  procedure  convergence,  and  the  final 
electric  field  distribution  E(r)  does  not  depend  on  the 
initial  approximation  of  E0(r).  However,  the  voltage 
drop  is  kept  constant. 


3.  Results 

It  was  found  that  the  iterative  process  is  converging 
only  in  the  case  of  moving  striations.  In  spherical 
discharge,  the  solution  for  ion  density  and  electric  field 
are  looked  for  in  the  form 
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where  n°j  and  E°  are  characteristic  values  of  ion  density 
and  electric  field  at  r=R,  ff.(/,r)  and  E(t,r)  are 

dimensionless  functions  presenting  waves  moving  with 
the  phase  velocity  vv/  = v^r//?)2^,  the  constant  v°j= 

otPiE0  has  the  order  of  ion  drift  velocity.  It  should  be 
noted  that  for  a<0.25  (standing  striations  corresponds  to 
a=0)  the  iterative  process  is  not  converged.  Below  all 
results  are  presented  for  a=0.5.  In  the  general  case, 
parameter  p can  be  in  the  range  P=l-2.  For  P=2,  average 
field  and  striation  velocity  vst  are  constant.  Ion  and 
electron  densities  vary  inversely  proportional  to  the 
square  of  the  distance  r from  the  anode.  Such 
dependencies  were  obtained  in  drift-diffusion 
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approximation  of  a spherical  glow  discharge  [3].  For 
P=l,  the  average  electric  Field  and  densities  ni5  nc  vary 
inversely  proportional  to  r.  In  this  case,  the  striation 
velocity  decreases  with  the  increase  of  the  distance  from 
the  anode. 

In  Fig.  1,  an  example  of  self-consistent  field  distribution 
and  electron  density  is  presented  for  argon  pressure 
p=0.5  Torr.  It  is  seen  that  average  electric  field  is 
modulated  by  several  non-equidistant  striation  peaks. 
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Fig.l.  Electric  field  and  electron  density  radial 
distributions  in  spherical  striated  glow  discharge  for 
p-1,  p=0.5  Torr. 


The  radii  of  striation  obey  the  geometric  series  law 
rn~ripn  (see,  Fig.  2).  Such  dependence  was  observed  in 
experiments  [1]  (however,  not  for  argon  but  for  the 
mixture  of  N2  with  small  addition  of  acetone). 


n 

Fig.2.  Dependence  of  striation  radii  rn  on  the  number  of 
striation,  n,  for  the  case  P=! . 


According  to  the  energy  balance  equation 

"T = - P-'O  ) - X Pr(r)  (6> 
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the  difference  between  energy  gain  from  the  electric 

field  Pr(tj  and  the  energy  loss  in  elastic  F*( r)  and  all 

inelastic  collisions  P^tfr)  is  compensated  for  by  the 

spatial  divergence  of  the  energy  current  density,  see 

Fig.3.  It  is  seen  that  there  is  a phase  shift  between  the 

maximum  of  energy  gain  in  the  field  and  maximum  of 


energy  losses  in  the  inelastic  collisions,  which,  in  fact, 
stimulates  the  propagation  of  waves  (moving  striations). 


Fig.3.  Normalized  energy  balance  terms.  ♦♦♦♦ 

divergence  term, gain  from  electric  field,  OOO 

loss  by  elastic  collisions,  -m-m-m-  loss  by  excitation 
E|=l  1 .27  eV,  a a a-  loss  by  excitation  E2=l  1 .7  eV,  solid 
line  - loss  by  ionization  Ej=  15.7  eV. 

4.  Conclusions 

A spherical  glow  discharge  has  been  studied  in  self- 
consistent  kinetic  consideration  of  moving  striations. 

It  is  shown  that  two  types  of  plasma  parameters 
distributions  are  possible  in  a stratified  spherical 
discharge:  1)  constant  averaged  field  in  the  spherical 
discharge  gap  and  electron  density  nc(r)~T*2  with  peaks 
due  to  striations,  and  2)  averaged  field  and  electron 
density  are  inversely  proportional  to  the  distance  from 
the  anode.  Both  regimes  were  observed  in  experiments. 
It  is  still  not  clear  why  real  system  in  PC  of  spherical 
glow  discharge  chooses  one  of  these  alternatives. 

The  presented  kinetic  model  is  applied  for  low-pressure 
inert  gases.  For  the  pressures  p>1.5  Torr  striation  are 
damped  due  to  energy  losses  in  elastic  collisions.  In 
these  cases,  it  is  necessary  to  take  into  account  the 
excitation  of  molecule  metastable  levels  and  radiative 
processes. 

Special  attention  should  be  paid  to  ionization  and 
recombination  processes,  which  equalize  creation  and 
loss  of  charged  particle  on  the  striation  length. 
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